JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 65B, No. 3, July-September 1961 
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King's integral expression governing the problem of the time-harmonic motion of a 
piston in a rigid wall is simplified. 



In a previous paper [4] 2 the case of an arbitrarily 
(not necessarily time-harmonic) moving disk in a 
rigid wall was discussed. The starting point was 
the eq (12) in the above quoted paper, 

(1) 

which, when y is replaced by ik (k = 2ir/\) repre- 
sents the velocity potential for the time harmonic 
case (King's integral [2]). The investigation of 
transient phenomena as carried out in [4] did not 
involve (1) directly, but its inverse Laplace trans- 
form with respect to y, and therefore not much 
attention was paid to (1) itself. But the computa- 
tion of the acoustic pressure distribution for the 
time-harmonic case makes an investigation of (1) 
itself imperative. Numerous approximate evalua- 
tions of (1) (with 7 replaced by ik, &>0) have 
been carried out. (for instance, [1], [2], [3], [5], [6]). 
It may be noted in this connection that (1) can be 
considerably simplified in as much as that it can be 
transformed into an integral expression with an 
elementary integrand which is extended over a 
finite interval of integration. However, the ex- 
pressions (22) and (23) in [4] for the transient field 
can be obtained as easily from this simplified ex- 
pression for U. 

We proceed now to replace the factor (t 2 +7 2 )~* exp 
[— 2(t 2 +y 2 )*] by the integral formula in [4(9)]: 

<T 2 + 7 2 )-*exp [-£(t 2 +t 2 )*] 

°° J (Tt)t(t 2 +z 2 )-i exp [-y(t 2 +z 2 )*] dt. 

Upon interchanging the order of integration, one 
obtains instead of (1) 

U=a i " t(t 2 +z 2 )~* exp [-y(t 2 +z 2 )*] 

i Jo(rp)Jo(Tt)Ji(Ta)dT \dt. 

The inner integral is known [4(20)] and therefore, 
instead of (1), 
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r* co 

U=a\ t(t*+z*)-icxp[-y(t 2 +z*)i]f(t)dt, (2) 



where 



/(OH 



(wa)- 



arccos 



t 2 +p 2 -a 2 
2pt ' 



0<a<\t-p\ 
\t-p\<a<t+ P 
a>t+p 



Again, as before, two cases A and B have to be 
distinguished. The results are 



Case A f a^>p: 

r a~ l 
/(*) = -\ Orft) -1 arccos 

Case B, a<Cp: 
f(t) = (Tra)~ 1 arccos 



t 2 +p 2 - 



2 P t 



t 2 +p 2 -a 2 
2pt 



0<t<a- P 

a— p^t^a-^r p 

t>a+ P 

p—a<^t<^p-^a 



=0 otherwise. 
Therefore, inserting/^) into (2), 
Case A, p<a; 



U-- 



P t(t 2 +z 2 )~ l/2 exp[-y(t 2 +z 2 ) 1/2 ]dt 

na+p 

+7T" 1 t{t 2 + Z 2 )- 1 ' 2 

J a-p 

exp[- 7 (Z 2 +2 2 ) 1/2 ] arccos C^f^ ) 
or, upon evaluating the first integral, 

U=- exp(-yz)-- exp{-y[z 2 +(a-p) 2 ] l/2 } 
7 7 



dt. 



+7T- 



t(t 2 +Z 2 ) 
J a-P 



-1/2 



exp[- 7 (< 2 +s 2 ) 1/2 ] arccos (^fy^V- 
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> s ( t2+ o p L a * ) dt - 



Case B, p^>a: 

J P-a 

! +2 2 ) 1/2 ] arccos . 

Substitute t 2 -\-z 2 =a 2 and get 

?7=iexp(-T^)-iexp{- 7 [^ 2 + (a-p) 2 ] 1/2 }+T- i r J 
7 7 J R' 

2 ( Q ,2_^2)l/2 J^ ^(^ ( 3 ) 

?7=7r _1 exp(— 70:) 

[-^2 _ 2 2 _|_ ^2 _ a 2-| 

Lw=W^J^' »>«■ (4) 



arccos 



As in [4(21)], 7? and R' denote the smallest and the 
largest distance of a point of observation from the 
circumference of the disk respectively. These prop- 
erties are given by 



/? = [z 2 + (a+ P ) 2 ] 1/2 , B' = [z*+{a- P n i i\ 

The expressions (3) and (4) for U are considerably 
simpler than [4(12)]. For the time-harmonic case, 
7 has to be replaced by ik. The expressions [4(22)] 
and [4(23)] for a "Dirac pulse" are immediately 
obtained from (3) and (4) by taking their Laplace 
inversions with respect to 7. 
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